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ABSTRACT 

We investigate the dust extinction properties in the host galaxy of GRB 050904 at z=6.29 by 
analyzing simultaneous broad band observations of the optical and UV afterglow at three different 
epochs. We show that the peculiar afterglow spectral energy distribution (SED) observed at 0.5 days 
and at 1 day after the burst (1.6 and 3 hours rest frame) cannot be explained with dust reddening 
with any of the extinction curves observed at low redshift. Yet, the extinction curve recently inferred 
for the most distant BAL QSO at z=6.2 nicely reproduces the SED of GRB 050904 at both epochs. 
Our result provides an additional, independent indication that the properties of dust evolve beyond 
z^6. We discuss the implications of this finding within the context of the dust production mechanisms 
through the cosmic ages. 

Subject headings: gamma-ray: bursts — dust, extinction — galaxies: high-redshift — galaxies: ISM 



1. INTRODUCTION 

Until recently, the investigation of the ISM at high 
redshift has been mostly based on systems observed in 
absorption along the line of sight of bright quasars or in 
Lyman break galaxies. However, the bright emission of 
Gamma Ray Bursts (GRBs) and the broad redshift dis- 
tribution of these objects (from local to z=6.3, so far) 
has highlighted GRBs as a new, independent tool to in- 
vestigate the ISM at high redshift. 

GRB 050904 was a bright, long burst occurred at red- 
shift z=6.29, the most distant GRB identified so far. 
Previous studies have shown that the progenitor was a 
massive sta r embedded in a dense, metal enriched molec- 
ular cloud (jCampana et al.ll2006t TFrail et al.ll2006f ). This 
burst presents two peculiarities with respect to lower red- 
shift GRBs. First, X-ray data at early times (within the 
first few minutes) shows a large column of gas along the 
line of sight, but little associated optical dust extinction 
(jCampana et al.l[2006| ). Note also that the gas absorption 
is found to decrease rapidly within the first few hours. 
Second, its afterglow presents a peculiar flux suppression 
at A re st ~ 1300 A (in addition to the Lya blanketing ef- 
fect), at 0.5 days after the burst, while the flux redward 
1600 A shows negligible extinction (lHaislip et al.|[2006h . 

In this work we further investigate these issues, and 
in particular the latter, in terms of dust properties. We 
collect broad band afterglow observations (from near-IR 
to X-rays) at three epochs after the burst, and fit the 
data with an intrinsic power-law and dust reddening. In 
addition to the 'standard' extinction curves typically as- 
sumed to model dust absorpti on in extragalactic objects , 
and in lower redshift GRB's (|Stratta et al.ll2004 120051 : 
iKann et a l. 2006), we also tested the extinction curve 
inferred from the most distant Broad Abso rption Line 
(BAL) QSO, at z=6.2 (jMaiolino et al.ll200l . We show 
that the latter extinction curve is in excellent agreement 
with the data of GRB 050904, providing further, inde- 
pendent evidence for an evolution of the dust properties 
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at z>6. We discuss the implications of our results also 
for what concerns the origin of dust at high redshift. 

2. MULTIBAND SPECTRAL ANALYSIS 

2.1. X-rays and intrinsic SED 

The X-ray spectra have been extracted in the 0.3-10.0 
keV energy range from Swift X-ray Telescope (XRT) data 
following standard procedures. At ~ 0.5 days after the 
burst trigger t (2005, Sept. 5.07785 UT), the X-ray 
emission shows a strong fl aring activity (jCusumano et a\\ 
l2007tlGendre et al.l[2006h . In order to avoid any possible 
spectral contamination from the flares in our continuum 
analysis, we did not take into account data in tempo- 
ral intervals where flares were present. We fit the 0.3- 
10.0 keV data by assuming an absorbed power law. We 
find evidence for little absorption in excess to the Galac- 
tic v alue (N§ al = 5 x 10 20 cm -2 , iDickev fc Lockmanl 



1990T). in agreement w ith the analysis at late times by 
Campana et al.l (|2006[ ). The best fit spectral energy in- 
dex is a = 1.2 ± 0.2 (F v cx v~ a ). This is significantly 
steeper than the spectral index measured just after the 
GRB event, suggesting that the synchrotron spectral 
break shifted toward lower energies between the first 
XRT observation, ~ 3 minutes after the GRB event, 
and the observations at 12 hours (in agreement with 
ICusumano et al.|[2007T ). 

The X-ray emission after 0.5 days decreases rather 
sharply, preventing a detailed spectral analysis after this 
epoch. Thus, we are able to directly compare the op- 
tical fluxes with the X-ray spectrum only at 0.5 days 
after the burst trigger. However, the spectral and tem- 
poral properties of this afterglow indicate that the cool- 
ing frequency is redw ard of the optical range already few 
hours after the burst (IKann et al.lf2007l : iFrail et alJl2006b 
ICusumano et al.ll2007[ K As a consequence, no spectral 
break is expected between the X-rays and the optical en- 
ergy range at any of the epochs of our analysis, nor the 
power-law index is expected to change in such a time in- 
terval. Therefore, we can safely assume that the intrinsic 
GRB spectrum at 0.5 days (observer frame) is a single 
power-law extending from the X-rays to the optical. At 
1 and 3 days we assume that the intrinsic SED has the 
same power-law index as at 0.5 days. 
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2.2. Optical-NIR 

In order to build a broad band SED and investigate 
its possible temporal variations, we collected multiband 
optical and near-IR photometry at three epochs: i +0.5, 
1 and 3 days after the burst trigger to (observer frame). 
These specific epochs were chosen because both near- 
IR and Z-band photometry (which are crucial to investi- 
gate the rest-frame UV SED) were obtained around these 
times. More specifically, at ^0.5 days after the burst we 
used Z, J, H photometry obtain ed with UKIRT an d K' 
photometry obtained with IRTF (jHaislip et al.ll2006l ). At 
~1 day after the burst we used J, H, Ks, and z band pho- 
tomet ry obtained with ESO/VLT by Tagli aferri et al.l 

(2005) . At ~3 days after the burst we used J, H, Ks 
photometry from lTagliaferri et al.l (120051) a nd z' obtained 
with Gemini South by IHaislip et all (120061) . 

Small temporal differences between the times of each 
individual observation for the three selected epochs were 
accounted for by assuming that the optical-UV after- 
glow declines with time as a broken power law. More 
specificall y, the best fit indices to the light curve ob- 
tained bv lTagliaferri et alJ (|2005f ) are 5 1 = -0.72 ± 0.17, 
&2 = —2.4 ± 0.4 a nd a break at tb = 2.6 ± 1 days after 
the burst (see also iKann et al.l (|2007f )). 

We corrected the observed magnitudes for the Galactic 
extinction toward the direction of this burst (E(B-V)= 
0.066) by assuming the lCardelli eTall (fl989h extinction 
curve with R v = A v / E{B -V)= 3.1. 

Fluxes in the z band were corrected for IGM Lya blan- 
keting effect by renorma l izing the afterglow spectrum 
obtained by iKawai et alJ ((2006) (at 3 days) so that its 
convolution with the specific z filter profile (including 
the detector response) yields the photometry observed 
at each epoch. From the re-normalized spectrum we 
derived the continuum flux density at 9300A (i.e. just 
redward of the damped wing of the IGM Lya absorp- 
tion), corresponding to A res t = 1275A. As we shall see 
in Sect. 13. li there is only little and marginal evidence 
for any variation of the (rest - frame ) UV SED between 
the epoch of the IKawai et a l. ( 2006) spectrum (3 days) 
and the previous epochs. This justfies the extension of 
the IGM Lya blanketing correction obtained at 3 days 
to the other epochs. Note that the z-band photometry 
in the three epochs is obtained with quite different filters 
at different telesc opes. In particular, the z filter (FORS2 
at VLT) used by iTagliaferri et alJ (|2005f ) at 1 day after 
the burst is much r edder than the z-band filters used by 
IHaislip et al.l (|2006h at the other two epochs. As a conse- 
quence, the z-band correction factors for Lya blanketing 
effect are significantly different at the three epochs, and 
specifically: a factor of 3.02 at 0.5 days, 1.27 at 1 day 
(i.e. nearly unaffected by IGM absorption), and 2.20 at 
3 days. 

3. DUST PROPERTIES 

3.1. Evidence for dust evolution at z>6 

At 0.5 days after the burst, a peculiar flux suppres- 
sion in the z b and ( X r e. K t ~ 1275A) has been detected by 
IHaislip et al.l (|2006| ). after accounting for IGM Lya ab- 
sorption. We find the same z band flux suppression also 
at 1 day after the burst, with a totally inde pendent set of 
data ( Fig. [J) . The flux drop is so sharp that Ha islip et all 

(2006) ruled out dust reddening as a possible cause. 



TABLE 1 

Best fit parameters of the GRB 050904 SED at different 

EPOCHS. 



Extinction 


A(3000A) a 


x"/d.o.f. b 


curve 


mag 




t — to = 0.5 day 


Galactic 


0.53 ±0.08 


38.7/(12-2) 


SMC 


0.33 ±0.06 


33.3/(12-2) 


Calzetti 


0.43 ± 0.06 


31.4/(12-2) 


QSO@z=6 


0.89 ±0.16 


15.1/(12-2) 


t — to = 1.0 day 


Galactic 


0.43 ±0.18 


27.9/(4-2) 


SMC 


0.48 ± 0.08 


17.4/(4-2) 


Calzetti 


0.70 ± 0.06 


24.3/(4-2) 


QSO@z=6 


1.33 ±0.29 


0.01/(4-2) 


t — io=3.0 day 


Galactic 


< 0.3(2cr) 


3.2/(4-2) 


SMC 


0.27 ±0.07 


1.1/(4-2) 


Calzetti 


0.40 ± 0.07 


2.2/(4-2) 


QSO@z=6 


0.46 ± 0.28 


0.6/(4-2) 



Wc further investigate the dust extinction scenario to 
tackle the issue of the peculiar z band flux suppression. 
Beside the standard extinction curves inferred in the lo- 
cal univers e (Galactic, SMC and starburst attenuation 
curve, from lCardelli et al.lll98l |Pelfl99l iCalzetti et all 
|1994[) 3 . we also consider the extinction curves computed 
bv IMaiolino et all (|2001l ) based on a dust model skewed 
toward large grains , which apply to some low-z GRBs 
(IStratta et al.l 12004 |2005|) . and the extinction curve m- 
ferred bv IMaiolino et al. |2004) for the most distant 
BAL QSO at z=6.19 (hereafter "QSO@z=6" extinction 
curve) , illustrated in Fig. O 

From the optical-to-X-ray spectral analysis, we find 
that at io+0-5 days the optical data require additional 
extinction in addition to the Galactic value to reproduce 
the flux values expected from a power-law with the same 
spectral index of the X-ray spectrum (see previous sec- 
tion). No 'standard' extinction curve can reproduce the 
observed SED with the peculiar flux suppression in the 
z band, as shown in the first three panels on the top of 
Fig. HJ The extinction curves from a dust model skewed 
toward larger grains also fail to reproduce the flux in 
the z band, although they provide a better fit to the J,H 
and K photometry than the 'standard' extinction curves. 
The only extinction curve that provides an acceptable fit 
is the QSO@z=6 one with A(3000A) - 1 mag (Fig. [Q 
rightmost top panel). The same result is obtained by 
using the independent set of data at to+1 day (Fig. [1] 
second row of panels). 

The statistics of the fits with the different extinction 
curves (along with the inferred absolute extinctions) are 
given in Table Q] It is clear that the improvement of the 
fit with the QSO@z=6 extinction curve is highly signifi- 
cant both at 0.5 and 1 days after the burst. 

This result is still true at 3 days, though with a lower 
significance. We also note that at the latter epoch there 
is marginal evidence for a decrease of extinction (Tab. [I}. 
This could be regarded as an indication of dust destruc- 
tion by the blast wave three days after the burst, or al- 
ternatively that at this late epoch the fireball has a pro- 



3 In the ca s e of the starburst attenuation curve we use the 
ICalzetti et~aTI l(l99lh expression for the UV continuum that, at 
3000A, is 0.4 times the extinction for the ionized gas. 
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Fig. 1. — optical-UV rest-frame SED of GRB 050904 at 0.5 (top), 1 (middle) and 3 (bottom) days after the burst (observer frame). The 
corresponding observed bands are marked on the top-leftmost panel. The source of data at each epoch are reported on the righ hand side 
(H06 = Haislip et al. 2006, T05 = Tagliaferri et al. 2005). At each epoch we show the best fit with each the four extinction curves shown 
in Fig. [2] as indicated on top of of each column. 



jected size larger than the obscuring cloud and therefore 
the emitting region is less absorbed. However, since the 
decrease of extinction at 3 days is only marginally signif- 
icant, we do not discuss this issue further and focus only 
on the extinction curve properties. 

Our result indicates that two totally different classes of 
objects at z > 6, QSOs and GRBs, are characterized by 
the same extinction curve, which is different from that 
observed at lower redshi ft. More specif i cally, the QSO 
at z=6.2 investigated bv iMaiolino et afl (|2004l ) and the 
GRB at redshift 6.3 investigated in this paper, provide 
independent evidences for an evolution of the dust prop- 
erties beyond z ~ 6. 

3.2. The nature of dust at z>6 

A transition in the properties of the dust at z ~ 5 — 6 
is theoretically expected. Indeed, the major source of 
dust in the local universe are the envelopes of AGB 
stars, which requ i re about 1 Gyr to evolve in large num- 
bers (lpwekll2005t iMorgan k Edmundsll2003l : iMarchenkol 
120061: iTodini k Ferraral |2001[) . At z>5 the age of the 
Universe is less than 1 Gyr, therefore AGB stars fall 



short of time to produce enough dust. As a con- 
sequence, galaxies at z>5 should be devoid of most 
of the dust, since the major source of dust (AGB) 
is missing. Nonetheless, significant masses of dust 
are still observed i n distant QSOs up t o z=6.4 (e.g. 
Bertoldi et ail 120031 : iPriddev et ail 120031: iRobson et all 

Theoretical models have 



2004 iBeelen et alJ [200l 



shown that an alternative source of dust on short time 
scales are core-collapse supernovae, which could there- 
fore be the prim ary dust production mechanism in 
the e arly un iverse (ITodini fc Ferrardl200lUNozawa et al 
I2003t lDwekll200a IMorgan fc Edmund s 2003; Marchenkc 



|2006[ ). From the observational point of view it is still not 
clear what is the actual efficiency of dust production in 
SN ejecta (ISugerman et all l2006t lErcolano et all 120071: 
Krause et all 12004 IWilson k, Batrlal 120051 : I Green et all 
20041) . However, both |Maiohno_eialJ (|2004( l 



and 



Hirashita et alJ (|2005l ) found a good agreement between 



the extinction curve expected from SN dust and the ex- 
tinction curve observed in the most distant BAL QSO 
at z=6.2 (QSO@z=6), suggesting that in the latter ob- 
ject dust is predominantly produced by SNe. The re- 
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Fig. 2. — Extinction curves tested in our analysis, normalized to 
A r cst = 3000A. The photometric bands are shown shifted to the 
corresponding central wavelengths at the rest frame of GRB 050904 
at z=6.29. The extin ction curve inferred for the QSO at z=6.2 
lIMaiolino et al.l I2004T) is directly observed at 1260A < A rost < 
3300A, while it is extrapolated beyond these limits by using the 
theoretical SN dust extinction curve. However the latter extrapola- 
tion is unimportant since for GRB 050904 at z=6.29 we essentially 
use the same observed bands as for the QSO at z=6.2. 

suit obtained in this paper indicate that also in the host 
of GRB 050904 dust is likely produced mostly by SNe. 
Altogether, these findings suggest that, more generally, 
regardless of the specific class of object, at z>6 SNe are 
the main source of dust. 

An independent evidence for different properties of the 
dust in the host of GRB 050904 comes from the compar- 
ison between the large column of absorbing gas observed 



at early times (Ng ~ 8 x 10 22 cm -2 , ICusumano et al.1 
l2007t [Campana et al.ll2Q06h and the little associated ex- 
tinction observed in the optical-UV rest-frame. More 
specifically, the latter implies an Ay /N# ratio more than 
50 times lower than the Galact ic value, which is unprece - 
dented in lower redshif t GRBs (jStratta et al.ll2004l2005h . 
ICampana et al.1 (|2006f) ascribe this effect to dust mostly 
composed of silicate grains (possibly produced by Pair In- 
stability SNe) which are destroyed soon after the burst. 
Here we note that, independently of the dust composi- 
tion, a strongly reduced dust-to-gas ratio is naturally ex- 
pected at z~6 as a consequence of the dust evolutionary 
time scales. Indeed, as discussed above, the lack of the 
important contribution to dust production from AGB 
stars at z>6, makes the dust-to-gas ratio in the early 
universe necessarily much lower than observed locally. 

Furthermore, we note that the variation of dust con- 
tent as a function of redshift may be at the origin of 
the shortage of optical GRB detections at lower redshift. 
Indeed, only ^50% of the Swift GRBs have shown an op- 
tical counterpart, even b y taking advantag e of the quick 
Swift GRB localization (|Fiore et al.l I2007T ). The larger 
dust content at z<6, expected by the dust evolution- 
ary scenarios, may prevent the optical detection of low-z 
GRBs. This issue will be discussed more extensively in 
a forthcoming paper. 
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